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Receptivity in parallel flows: an adjoint approach 

By D. C. Hill 


1. Motivation and objectives 

Linear receptivity studies in parallel flows are aimed at understanding how ex- 
ternal forcing couples to the natural unstable motions which a flow can support. 
The vibrating ribbon problem (Gaster (1965)) models the original Schubauer and 
Skramstad (1947) boundary layer experiment and represents the classic bound- 
ary layer receptivity problem. The process by which disturbances are initiated in 
convectively-unstable jets and shear layers has also received attention (Balsa (1988), 
Huerre and Monkewitz (1985), Tam (1978)). 

Gaster (1965) was the first to handle the boundary layer analysis with the recog- 
nition that spatial modes, rather than temporal modes, were relevant when studying 
convectively-unstable flows that axe driven by a time-harmonic source. The ampli- 
tude of the least stable spatial mode, far downstream of the source, is related to 
the source strength by a coupling coefficient. The determination of this coefficient 
is at the heart of this type of linear receptivity study. 

Traditionally, the Briggs method is applied to Fourier- inversion integrals to find 
the asymptotic temporal and spatial behavior after a time-harmonic source has 
been switched on. Ashpis and Reshotko (1990) give a detailed description of this 
procedure for the vibrating ribbon problem. Unfortunately, the coupling coefficient, 
which gives a measure of the amplitude of the asymptotic response relative to the 
amplitude of the source, does not have a very convenient form, either from the point 
of view of interpretation or of computation. The expression involves derivatives 
with respect to wavenumber of either a flow quantity or the dispersion relation 
(Ashpis and Reshotko (1990), Kozlov and Ryzhov (1990), Balsa (1988), Huerre and 
Monkewitz (1985), Gaster (1965)). 

Earlier work (Hill (1992)) indicates that adjoint eigensolutions for the global 
temporal modes of a cylinder wake characterize rather simply how receptive the 
wake is to control forces. Can this approach be adapted for the spatial instabilities 
that occur in boundary layers, jets, and shear layers? The first objective of the 
present study was to determine whether the various wavenumber derivative factors, 
appearing in the coupling coefficients for linear receptivity problems, could be re- 
expressed in a simpler form involving adjoint eigensolutions. Secondly, it was hoped 
that the general nature of this simplification could be shown; indeed, a rather 
elegant characterization of the receptivity properties of spatial instabilities does 
emerge. The analysis is quite distinct from the usual Fourier-inversion procedures, 
although a detailed knowledge of the spectrum of the Orr-Sommerfeld equation is 
still required. Since the cylinder wake analysis proved very useful in addressing 
control considerations, the final objective was to provide a foundation upon which 
boundary layer control theory may be developed. 
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2. Background 

2.1 The Lagrange identity 

The cornerstone of this work is the Lagrange identity (Ince (1944)), from which 
the adjoint equations are extracted. This remarkably powerful relation is implicit 
in the work of Salwen and Grosch (1981) (here onwards referred to as SG), and Sal- 
wen (1979), although it is never stated explicitly. As has already been intimated, 
the origins of the adjoint analysis can be traced back to work of Lagrange in the 
late 18th century (See Lagrange’s collected works (1867)). Although Fuchs (1858) 
wrote in German, Ince (1944) credits the use of the word ’adjoint’ to him. Adjoint 
eigensolutions play a key role in the imposition of solvability conditions (the Fred- 
holm alternative) and as such are indispensable to researchers studying nonlinear 
phenomena using bifurcation theory (Iooss and Joseph (1980)). 

The Lagrange identity is developed with reference to the linearized Navier- Stokes 
equations. If F(r) is a steady incompressible viscous flow field, then linear velocity 
disturbances v(r y t) and pressure disturbances p(r, t) upon this flow satisfy 


^ + L(V-,R)v + Vp = 0, 

V.u = 0, 

where the ith component of the linear operator L{V_\ R ) is 

/r/T _ \ <9 Vi 1 d 2 V t 


( 1 ) 

( 2 ) 


(3) 


For any pair of fields (v,p) and (v,p), (( v,p ) does not have to satisfy (1) and (2)), 
defined over the flow domain, the following Lagrange identity is satisfied: 


Fl'i) 

+L(V-,R)v + Vp).v + V.vp 


+ 


dv *v 

v. (-7^ + L{Vj, R) v + Vp) + pV.v 


= + V.J((v,p),(v,p)), 


(4) 


where L(V_\ R) is the adjoint linearized Navier-Stokes operator with components 


The vector J((u,p), (v,p)) is the bilinear concomitant with components 

(l{(v,p),(v,p))^ = Via,j + (T t jVi, 


(5) 


( 6 ) 


Receptivity in parallel flows 


229 


where 


r 1 <k>i , „ 

~ pSt} Rd Xj +VjV " 

( 7 ) 

~C 1 S5i 

= pS ‘ i+ Rd,/ 

( 8 ) 


Examining the second term in square brackets on the left hand side of the La- 
grange identity (4), we define the adjoint equations 

-£+t(V;R)v + Vp = 0, (9) 

V.5 = 0. (10) 

2.2 Bi- orthogonality 

A brief summary will be given in this section of the procedure by which bi- 
orthogonal spatial-eigensolution sets can be constructed for the Orr-Sommerfeld 
equation. This is essentially a review of the work of SG. 

Schensted (1960), and Drazin and Reid (1981) show how general disturbances 
in spanwise-bounded plane parallel flows can be expanded as a sum of temporal 
eigensolutions of the Orr-Sommerfeld equation. These eigensolutions together with 
a set of eigensolutions of the corresponding adjoint equation form a bi-orthogonal 
set, under the action of an appropriate inner product. SG develop this further 
and demonstrate how expansion of arbitrary disturbances as a sum of temporal 
or spatial eigensolutions may be carried out in unbounded flows where, in fact, a 
continuum of eigensolutions forms part of the spectrum. 

Following SG, let V (r) = U(y)x define a two-dimensional parallel flow in the 
xy-plane, for y x < y < y 2 . For a boundary layer we may have yi = 0 and y 2 = oo, 
whereas for a shear layer y x = -oo, and y 2 = oo. Writing v = V x (t/>(x, y, t)z ) for 
some stream function xj)(x,y,i), z being the unit vector normal to the plane of the 
flow, the governing equation is given by the z-component of the curl of (1) 




dy 2 dx R V 

Similarly, writing £ = V x (ij>(x,y,t)z), the adjoint equation is 

/ d „ d \ , - dU d 2 ^ 1 . - 

(dt + U dx) V * + 2 dy d-rSy + R V ^ = °' 

No-slip boundary conditions 


(ii) 


( 12 ) 


, dxb - dxh 

ip = — = 0 t and 0 = — = o 


( 13 ) 
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are imposed on any walls. As |y| - oo, disturbances are required to decay to zero. 

As a direct consequence, the {/-component of J approaches zero as y approaches the 
limits of the flow domain. 

With a, a being wavenumbers and being frequencies, let 

</>(*, y,i) = Ku{yy {wt+ax \ t A(x,y,t) = ^w(y)e-’ (i " +5l) (14) - 

be solutions of equation (11) and the adjoint (12), respectively (i.e. 4> au (y) satisfies ■ 

the Orr-Sommerfeld equation, and hu,(v) an adjoint Orr- Sommerfeld equation). 

An auxiliary problem can then be solved to obtain the pressure eigensolution and 
the adjoint 

p(x,y,t) = Paw(y) e * (u;t+0fX \ p( x ,y>t) — Paw(y) e * + (!5) 

Substituting these velocity and pressure fields into the Lagrange identity (with 

V(r) = U{y)x ), the entire left hand side is zero, from which it follows (SG, Salwen ~ 

(1979)) that _ _ ^ 

(u> — fe’) + (o — d)[<^ Qu; , <t>aus] = 0) (1®) " 

with 

<*»,*«> = ^ + ( 17) 

/ {(^a + 

Jy\ K 

/ ^ « j d 2 4>au aw \ _ 

+ U + oa + d )<pawdaw + raw ^~~dy <9y / 

+ + (18) = 

When considering the spatial stability problem, w is chosen to be purely real, 
and from the Orr-Sommerfeld equation, a set of characteristic values for a can then 
be found. This set can consist of discrete and continuum modes. For the adjoint 
problem, likewise selecting u = ui, a set of possible values for a may be obtained. 

From (16), we have 

(a - d)[</> aw , 4>aw] = 0, (1®) 

from which it follows that 

if [i«, ^ 0 , then a = a. (20) 

With appropriate care taken in handling the continuum modes, SG proceed to con- 
struct a bi-orthogonal eigensolution set in which for every solution <f> aw (y)e t(u " +ax) 
to (11), there corresponds a normalized adjoint solution <£aw (y) e ^ ^ = 

such that, for example, for discrete modes with wavenumbers a,', 

[$a,-w , ^aj w] &ij' 


( 21 ) 
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A similar result is quoted by SG for the case of temporal modes, based on the inner 
product (*,*). 

For an arbitrary solution s = (y(x, y),p(x, y))e twi of the linearized equations, the 
adjoint solutions can be used to decompose the disturbance field. For a chosen 
eigenmode, if 3 = (v au; (y),Pau;(y))e _, ^ wt+Qfr ^ is the adjoint eigenmode, the filtering 
operation 

f £-l{§.,l)dy (22) 

Jy\ 

returns the modal amplitude. 

3. Accomplishments 

In this section I will give examples of how the adjoint problem provides a refor- 
mulation of various aspects of the receptivity problem. 

3.1 Failure of hi- orthogonality 
It can be demonstrated that 


duj 

da 


[^cruM ^au>] 
{^orwi ^aw) 


(23) 


which I believe to be a previously- unknown form for the group velocity of Orr- 
Sommerfeld eigenmodes. This form has implications for the formation of orthogo- 
nad solution sets. In the event that the chosen eigensolution is at a branch point 
singularity, for which du)/da = 0, then clearly 




(24) 


(the denominator in (23) is always bounded). The normalization condition (21) 
would require that the bilinear operation (24) take the value 1. The failure of (21) 
is connected to degeneracy of the wavenumber eigenvalues as the branch point is 
approached in the complex-u> plane. 

3.2 The vibrating ribbon problem 

The classical vibrating ribbon problem was discussed in some detail recently by 
Ashpis and Reshotko (1990). A vibrating ribbon is placed at the wall beneath a 
parallel boundary layer flow, and, after t = 0, it oscillates with velocity 

V h (x,t) = 6(x)e ,ut y. (25) 


The classical result indicates that a long time after the initiation of the excitation 
and far downstream of the ribbon, the y-ve locity of disturbances is 


v{x,y, t) 


o(jr,a,v) i (ut +*s) 


(26) 
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where a(u>) is wavenumber of the fastest growing discrete spatial mode of the Orr- 
Sommerfeld equation at frequency w. The field t>(y;a,u>) is the y- velocity field of 
the mode. The factor dv/da arises from contour integration around a pole in the 
complex wavenumber plane. 

Beginning with (26), by differentiating the linearized equations of motion with 
respect to a and making use of the Lagrange identity, I have shown that 

v(x,y,t)= p aul (0) v(y;a,u>)e l( - wt+at) . (27) 

Details of the analysis are not given here. The adjoint pressure at the wall, p a w(0), 
gives the amplitude of the response to the ribbon oscillation. 

8.8 Excitation of a free shear layer 

Huerre and Monkewitz (1985) consider the response of an inviscid incompressible 
free shear layer ( U(y ) defined for -co < y < oo) to excitation by a point vorticity 
source positioned at :r = 0, y = yo, 

n(x,y,t) — 6(x)6(y - y 0 )e ,ut . (28) 


Quoting their result, far downstream the stream function behaves as 


rp(x,y,t) = 


( yo)4> QUf (: V ) 


c t(w*+ax) 


(29) 


(al7(yo) + u>)f£(y 0 ;a,u>) 
where <£ auJ (y ) is the least stable spatial eigensolution of the Rayleigh equation, and 

d$- 


D(y 0 ;o,w) = ($ + - 


dy 


•-£) ■ 
oy ' y=yo 


(30) 


Here $+(y;a,u) and $~(y;a,u;) are the upper (y > y 0 ) and lower (y < y 0 ) eigen- 
solutions of the Rayleigh equation, which decay as |y| — ► oo. 

I have taken (29) and simplified it to obtain 

xl>(x,y,t)= ^aw(yo) ^(y)e ,(w<+ax) . (31) 


In this case, the modal amplitude is simply the value of the adjoint stream function, 
^aw(yo), at the source location. 

3.4 Response to a general harmonic source distribution 

Consider a source distribution oscillating harmonically at real frequency in- 
cluding sources of momenta of strength y(r;cj), and mass sources ^(r;u?). To keep 
the analysis as general as possible, we consider a boundary layer flow in the region 
y > 0 so that boundary velocities u^,(x;u?) can be specified at y = 0. The governing 
equations are then 


uov + L(V\ R)v + Vp = q(r ; u>), 

V.u = <p(r] lo), 


(32) 

(33) 
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with 

v = on y == 0. (34) 

The sources are assumed to be localized in so far as they disappear for |x| larger 
than some value, say X > 0. 

Let the discrete spatial mode of interest have stream function <£au;(y) e, ^ <+orI ^ 
Typically the mode with largest value of — /m(a(u;)) would be considered. We wish 
to determine the amplitude of this mode far downstream as a result of the excitation 
by the various sources. 

For a convectively-unstable mode which is traveling downstream, the stream func- 
tion behaves as 

a *«4y)«* M+< “ ) (35) 

in the region downstream of sill sources. The amplitude a is given explicitly by 


fx 2 r oo too 

a- j / g(r,u).v QU> (y)e~' az dydx + / / <p(r\u>)p au (y)e~ ,ax dydx 

J X \ J 0 J X\ JO 


+ 


/ 1>j(z;u>).5 0 

Jxx 


(36) 


,e- iax dx, 


where x\ < —X,x 2 > X. Clearly, if there are no sources, then a = 0. The field 


Uaw(y) = ( 37 ) 

is the adjoint velocity field, p QW (y) is the corresponding adjoint pressure, and the 
complex vector 

S ~"(*“* + 5T!f f ),.. (38) 

is the adjoint stress . 

The effect of some other types of sources can also be deduced. For sources of 
vorticity, the term involving £(r;u>) is rewritten as 


rX2 r OO _ rX2 y*0O 

/ / £(i;w).V x {xj)z)dydx = / / ft(r;o;)<£ OU ;(y)e “ ,az dydx (39) 

J z\ Jo Jx i Jo 

where f)(r;u>) = z.(V x £(r;u>)) is the vorticity source distribution which would ap- 
pear, for example, on the right-hand side of the Orr-Sommerfeld equation. Vorticity 
sources in the flow are weighted by the adjoint stream function. 

Rather than specify a velocity at the wall (i.e. at y = 0), suppose that the wall 
is in motion, oscillating about its mean position with a small velocity p rf (:r;u;). 
Linearizing this boundary condition, it follows that the boundary integral in (36) 
can be re-expressed as 
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where 



-(( 


Pau> 4 " 


i dU du aul \ „ 
u;i? dy 5y / 


1 dug* 

R dy /y=o 


(41) 


we will call the modified adjoint stress. 

In each instance, the streamwise integration is weighted by e~' ax . Since, typically, 
e iax g rows downstream, e~ iaz will grow upstream. There is no surprise here since 
sources further upstream will have a greater contribution to the far field disturbance 
amplitude; the response to such sources has convected further and hence has grown 


more. 


The deductions made in sections 3.2 and 3.3 can now be reconfirmed. For the 
vibrating ribbon problem, the response to boundary motion (25) (ri < 0, x 2 > 0) 
predicted by (36) is 


f 


6(x)y.S au e~ iax dx = p ou ,(0). 


(42) 


For excitation of a free shear layer by a vorticity source (28) with the y integration 
in (39) now extending from — oo to oo, the amplitude of the response is 


a: 


6(x)6(y - yo)4>aw(y)e ,az dydx = f> aLJ (yo)- 


(43) 


The amplitude of a particular spatial eigensolution generated by an upstream 
time-harmonic source distribution can be expressed as a weighted integral of the 
sources whether they are within the flow or upon a flow boundary. The weighting 
functions are simply different field quantities of the adjoint eigensolution corre- 
sponding to the mode being considered. The field quantity which is appropriate 
depends upon the nature of the source. The following table summarizes the cases 
considered in this section. 

Source type Adjoint weighting factor 


description 

symbol 

symbol 

description 

integrand 

momentum 

2 .(r;u>) 

£aw(y)e- ,ax 

velocity 

q(r-,u).v au {y)e- ,ax 

mass 

V?(r;u;) 

p a u l (y) e ~ ,az 

pressure 

<f(r\uj)p au/ (y)e-' ax 

vorticity 

Q(r; u) 

4>Uy)e- ,ax 

stream function 

ft(r;ui)<j>aw{y)e- ,az 

velocity at 
boundary 

v b (x-,u) 

S Q u e- iax 

adjoint stress 

v b (x\uj).S au ,e-' ox 

velocity of 
boundary 

v d (x;u>) 

£Le-' ax 

modified 
adjoint stress 

v d (x]u).Sl au! e-' oz 


For momentum, mass, and vorticity sources, the integration is made over the entire 
flow domain in which sources are present. The resulting value gives the amplitude 
of the mode far downstream. For boundary sources, the integration is made over 
the boundary. 
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4. Future plans 

The next step in this work will be to obtain simple numerical solutions of the 
adjoint fields in flows such as the Blasius boundary layer. A map of the receptivity 
characteristics for these flows can then be found. This will supplement analytical 
studies of boundary layer receptivity (Goldstein (1983), Goldstein et al. (1983)). 

The coupling of free stream disturbances to boundary layer motions as a con- 
sequence of surface roughness is an important receptivity path (Goldstein (1985)). 
This has not yet been considered in the present work, and efforts will be made to 
extend the analysis to handle this scenario. 

In the area of control, a means of analyzing boundary layer control strategies 
will be pursued. Suppression of a global (temporal) instability, such as occurs in a 
cylinder wake, can be achieved by a small permanent alteration of the flow field. The 
corresponding spatial problem is more complex since the control forces in practice 
are localized in space (for example, a region of suction in a boundary layer (Saric 
and Nayfeh (1977)), with their effect being felt both upstream and downstream. 

Although reduction of the spatial growth rate of instabilities may be important, 
consideration will also be given to the alteration in the receptivity characteristics 
as a consequence of the presence of the control system. This can be quantified by 
examining changes in the adjoint field as a result of the control. 

The global instability problem for strongly non-parallel flow has already been 
handled successfully (Hill (1992)). After the control of spatial instabilities in parallel 
flows has been fully investigated here, it will remain to consider the spatial problem 
in non-parallel flows. 

It would seem inevitable that a connection will be established with the work of 
Herbert and Bertolotti (1987) on the Parabolised Stability Equations. Their studies 
on the evolution of a disturbance amplitude in a slowly evolving non-parallel flow 
would appear to be intimately connected to the present work, though they do not 
consider the receptivity problem explicitly. 

In the longer term, a study will be carried out of the crossflow instability on an 
infinite swept airfoil leading edge. This is a phenomenon of major technological 
importance, and it is hoped that a systematic means of analyzing control possibili- 
ties may be found. Providing the ability to analyze how secondary instabilities and 
turbulent flows respond to control forces also remains a long term goal. 
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